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Abstract 

Background: Analyses of mitochondrial (mt) genome sequences in recent years challenge the current working 
hypothesis of Nematoda phylogeny proposed from morphology, ecology and nuclear small subunit rRNA gene 
sequences, and raise the need to sequence additional mt genomes for a broad range of nematode lineages. 

Results: We sequenced the complete mt genomes of three Ascaridia species (family Ascaridiidae) that infest chickens, 
pigeons and parrots, respectively. These three Ascaridia species have an identical arrangement of mt genes to each 
other but differ substantially from other nematodes. Phylogenetic analyses of the mt genome sequences of the 
Ascaridia species, together with 62 other nematode species, support the monophylies of seven high-level taxa of the 
phylum Nematoda: 1) the subclass Dorylaimia; 2) the orders Rhabditida, Trichinellida and Mermithida; 3) the suborder 
Rhabditina; and 4) the infraorders Spiruromorpha and Oxyuridomorpha. Analyses of mt genome sequences, however, 
reject the monophylies of the suborders Spirurina and Tylenchina, and the infraorders Rhabditomorpha, 
Panagrolaimomorpha and Tylenchomorpha. Monophyly of the infraorder Ascaridomorpha varies depending on the 
methods of phylogenetic analysis. The Ascaridomorpha was more closely related to the infraorders Rhabditomorpha 
and Diplogasteromorpha (suborder Rhabditina) than they were to the other two infraorders of the Spirurina: 
Oxyuridorpha and Spiruromorpha. The closer relationship among Ascaridomorpha, Rhabditomorpha and 
Diplogasteromorpha was also supported by a shared common pattern of mitochondrial gene arrangement. 

Conclusions: Analyses of mitochondrial genome sequences and gene arrangement has provided novel insights into 
the phylogenetic relationships among several major lineages of nematodes. Many lineages of nematodes, however, are 
underrepresented or not represented in these analyses. Expanding taxon sampling is necessary for future phylogenetic 
studies of nematodes with mt genome sequences. 
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Background 

Nematodes (also called roundworms) are an extremely 
diverse group of bilateral animals with an estimate of 1- 
10 million species although only -25,000 species have 
been described [1]. Nematode taxonomy has traditionally 
been formed from morphology and ecology, e.g. the early 



* Correspondence: rshao@usc.edu.au; xingquanzhu1@hotmail.com 
+ Equal contributors 

3 Genecology Research Centre, University of the Sunshine Coast, Queensland 
4558, Australia 

1 State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of 
Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research 
Institute, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu Province 
730046, People's Republic of China 

Full list of author information is available at the end of the article 

(3 BioMed Central 



systems proposed by Schneider (1866) [2] on somatic 
musculature and by Cobb (1919) [3] on stoma armature, 
and the system proposed by Filipjev (1929) [4] on the 
presence and absence of zooparasitism. Chitwood (1937) 
[5] proposed initially the bipartite phylogenetic and taxo- 
nomic system based on the presence and absence of phas- 
mids; this system was further elaborated later by Maggenti 
(1963, 1983) [6,7] based on pharyngeal structure and ex- 
cretory systems. However, due to the lack of homologous 
characters and informative fossil records, and the exten- 
sive convergent evolution, it has been extremely difficult 
to derive a consistent phylogenetic and taxonomic frame- 
work for the phylum Nematoda from morphological and 
ecological characters [8-12]. In the past two decades, 
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nematode phytogeny and taxonomy have been revised 
with analyses of the nuclear small subunit (SSU) rRNA 
gene sequences [12-14]. The current working hypothesis 
of Nematoda phytogeny incorporated evidence from 
morphology, ecology and SSU rRNA gene sequence ana- 
lyses [13,14]. 

The mitochondrial (mt) genomes of bilateral animals 
usually contain 37 genes and a control region on a circu- 
lar chromosome, -16 kb in size [15-17]. Sequences of 
individual mt genes and whole mt genomes are widely 
used to infer phylogenetic relationships among animals 
at different taxonomic levels [18,19]. In a number of 
cases, arrangement of mt genes has also been used to re- 
solve long-standing phylogenetic relationships that could 
not be resolved by other means [20-22] . 

Recently, mt genome sequences have also been 
analyzed to understand the phylogenetic relation- 
ships among nematodes. Kang et al. [23] inferred the 
Nematoda phytogeny with the mt genome sequences of 
25 species; Park et al. [24] and Sultana et al. [25] ex- 
panded further the analysis to include 36 and 41 species 
respectively. A major difference between mt genome 
phytogenies and the current working hypothesis of the 
Nematoda is on the relationship among four infraorders: 
Ascaridomorpha, Rhabditomorpha, Spiruromorpha and 
Oxyuridomorpha. In the current working hypothesis, 
Ascaridomorpha, Spiruromorpha and Oxyuridomorpha 
are in the suborder Spirurina whereas Rhabditomorpha 
is in the suborder Rhabditina [13,14]. Further- 
more, Ascaridomorpha is most closely related to 
Spiruromorpha, whereas Oxyuridomorpha is sister to 
Ascaridomorpha + Spiruromorpha + Rhigonematomorph. 
In the phytogenies inferred from mt genome sequences, 
however, Ascaridomorpha is most closely related to 
Rhabditomorpha; in most analyses, Oxyuridomorpha 
is sister to Ascaridomorpha + Rhabditomorpha, whereas 
Spiruromorpha is sister to the group that contains 
Ascaridomorpha, Rhabditomorpha and Oxyuridomorpha 
[23,24]. Due to limited taxon sampling, however, the novel 
phylogenetic relationships inferred from mt genome se- 
quences in these studies need to be interpreted with 
caution and further test with more taxa from a wide range 
of lineages is necessary. 

In this study, we sequenced the complete mt genomes 
of three Ascaridia species in the family Ascaridiidae of 
the infraorder Ascaridomorpha. Ascaridia species are 
among the most prevalent and pathogenic parasitic 
nematodes found in domestic and wild birds and have a 
worldwide distribution [26]. These Ascaridia species 
share an identical arrangement of mt genes to each 
other but differ substantially from those of other nema- 
todes. We inferred the phylogenetic relationships with 
the complete mt genome sequences of the Ascaridia 
species and 62 other nematode species that have been 



sequenced to date. Our analyses support the close rela- 
tionship between Ascaridomorpha and Rhabditomorpha, 
and provide novel insights into other phylogenetic rela- 
tionships in the Nematoda. 

Results 

Mitochondrial genomes of three Ascaridia species 

The complete mt genomes of A. galli, A. columbae and 
Ascaridia sp. (GHL-2012) were 13,977 bp, 13,931 bp 
and 13,862 bp long, respectively (GenBank accession 
numbers JX624728, JX624729 and JX624730). Each mt 
genome contains 12 protein-coding genes, 22 tRNA 
genes and two rRNA genes on a circular chromosome 
(Figure 1). As in most other nematodes, atp8 gene is not 
present in the mt genome of A. galli, A. columbae and 
Ascaridia sp. (GHL-2012); all mt genes are transcribed 
from the same direction. TTT is used as the initiation 
codon in nad6 gene in the mt genome of A. galli. TTT 
was previously reported to be the initiation codon for 
cox2, cytb, nad4, nad3, nad2 and cox3 genes in the mt 
genome of Strongyloides stercoralis (Strongyloididae), 
but not for any other Ascaridida nematodes. All of the 
12 protein-coding genes have complete termination 
codons, TAA or TAG. The tRNA genes in mt genomes 
of A. galli, A. columbae and Ascaridia sp. (GHL-2012) 
range from 51 to 76 bp. The tRNA-Ser (AGN) and tRNA- 
Ser (UCN) have truncated secondary structures; these two 
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Figure 1 The mitochondrial genomes of three Ascaridia species. 

All genes are on the same DNA strand and are transcribed clockwise. 

Protein-coding and rRNA genes are indicated with the standard 

nomenclature. tRNA genes are indicated with the one-letter code of 

their corresponding amino acids. There are two tRNA genes for 

leucine: U for codons CUN and L 2 for UUR; and two tRNA genes for 

serine: S] for codons AGN and S 2 for UCN. "NCRL" refers to the large 

non-coding region. "NCRS" refers to a small non-coding region. 
v> ) 
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tRNAs lack a DHU stem but possess a T^C loop. In all 
other 20 tRNA genes, the T^FC arm and variable loop 
are replaced with a TV replacement loop. There are two 
non-coding regions (NCR) in the mt genomes of A. galli, 
A. columbae and Ascaridia sp. (GHL-2012) (Table 1). 



The longer NCR (NCRL) is between trnC gene and trnN 
gene, and is 610 bp (A. galli), 563 bp (A. columbae) and 
566 bp [Ascaridia sp. (GHL-2012)] respectively. The 
shorter NCR (NCRS) is between nad4 and trnM gene, 
and is 157 bp (A. galli), 91 bp (A. columbae) and 101 bp 



Table 1 List of annotated mitochondrial genes and regions of Ascaridia galli, Ascaridia columbae and Ascaridia sp 



Gene/region 




Position/length (bp) 






Start/stop codon 


Anticodons 


A. galli 


A. columbae 


Ascandia sp. 


A. galli 


cnlumhnp 


AsrnriHin 


coxl 


1-1 563 (1 563) 

1 1 JVJJ \ \ JVJjy 


1-1 563 (1 563) 


1-1 581 (1 581) 


ATG^AA 


ATG^AA 


GTG^AG 




tRNA-Cvs (C) 

ii\iN/ \ v_ y o V^—y 


1 565-1620 (56) 

I JVJJ I \J Z-\J \~J\JJ 


1 565-1620 (56) 


1 562-1 616 (55) 








GAT 


Nnn-rnHinn rpninn (NO 1 

1 NUI 1 V_ VJH 1 1 \KJ \ \Z\-jl\J\ 1 \ 1 N v — 1 — y 


i uz. i Z.Z.JW \vj i vj y 


1623-2183 (563) 

I UZ.J Z_ I D~J \~J\J^JJ 


1617-2182 (566) 

I \J I / Z_ I KJZ- \~J VJVJ ) 










tRNA-Asn (N) 


2231 -2284 (54) 


2184-2241 (58) 

Z_ 1 U 1 Z_Z_ i 1 


2183-2239 (57) 








G^ 


tRNA-Tvr (Y) 


2286-2344 (59) 
zzou /.jtt 


2246-2303 (58) 

ZZ^vJ Z.JUJ W t - ) / 


2240-2295 (56) 








GTA 


nod 1 


2342-3217 (876) 

ZJ^Z _>Z I / \(J / KJ) 


2301-31 76 (876) 


2298-3167 (870) 

ZZVU J I vj/ V<_)/ 


TTG/TAA 

I I VJ/ I / \/\ 


TO^AA 


G^^AG 




Qtp6 


3217-3813 (597) 


31 76-3772 (597) 


31 70-3766 (597) 

J 1 / v -J / VJVJ \.J -/ / y 


ATA/TAA 


ATA^AA 


^G^AG 




tRNA-l vs (K) 

ii\iN/ \ i_ y j \'v 


3816-3876 (61) 

J (J I VJ JU / U \\J I J 


3790-3R53 (64) 


3772-3833 (62) 










tRNA-l pu UUR (1 ■,) 

i\\\\r\ lcu V-2/ 


3876-3931 (Zft) 


3853-3907 (55) 


3831-3886 (56) 

JUJ I JUOU \J\J) 








TAA 


tRM A-Ser AGN fSO 


3932-3982 (51) 


3908-3962 (55) 


3887-3939 (53) 








GCT 


nod2 


3986-4828 (843) 

~J J\J\J \KJZ-KJ J 


3g nn _4808 (843) 

-J J\J\J \D\JD \D\~J J 


3943-4779 (837) 


1 1 VJ/ 1 / \/\ 


TO^AA 


TO^AG 




tRNA-lle (1) 


4832-4891 (60) 


4823-4884 (62) 


4781-4841 (61) 








GAT 


tRNA-Arg (R) 


4896-4950 (55) 

1 VJ v VJ 1 2Ju \~) ~) } 


4888-4944 (57) 


4845-4900 (56) 








ACG 


tRNA-GIn (O) 

ii\iN/ \ vj 1 1 i 


4951-5005 (55) 

"-JJ 1 -J VJ VJ J \~J ~J J 


4945-4998 (54) 


4901-4954 (54) 










tRMA-Asn (D) 

U\IN/\ / \-J kJ \LV J 


5005-5063 (59) 

JwwJ J VJVJ J J } 


501 1-5070 (60) 


4971-5025 (55) 








GTC 


tRNA-Glu (E) 


5065-5125 (61) 


5069-5126 (58) 


5025-5083 (59) 










/7T?S 


5123-5823 (701) 

J I z_ J J(_)Z_J \f \J \ J 


51 24-5826 (703) 

~J \ Z_ i JUZ.U \/ \J ~J ) 


5084-5774 (691) 

.J VJkJ" ~J / / 1 \\J J I J 










tRNA-Spr UCN (Sn) 


5826-5880 (55) 


5837-5891 (55) 


5781-5835 (55) 








TGA 


tRNA-Phe (F) 


5883-5939 (57) 


5898-5953 (56) 


5841-5897 (57) 








GAA 




5964-7067 (1 1 04) 

J J\J\ / VJVJ l \ \ \ \J\j 


5975-7075 (1101) 


5919-7019 (1 101) 


ATG/TAA 

r \ \ vj/ 1 r\i\ 


G^/TAA 


G^^AG 




tRNA-l pu CUN fl 0 

ii\iN/\ Lcru V— 1/ 


7067-71 22 (56) 


7077-71 34 (58) 


7020-7075 (56) 

/ wZw / KJ / U \~)\JJ 








TAG 


cox3 


7144-7914 (771) 


71 54-7918 (765) 


7097-7864 (768) 


ITTG^AG 


TO^AA 


^G^AA 




tRNA-Thr (T) 

LI\IN/\ MM V'/ 


7889-7944 

/ VJCJ _y / J V \ \~J\JJ 


7899-7956 (58) 


7845-7900 (56) 








TAG 


nod4 


7939-91 74 (1 236) 


7951-9186 (1 236) 


7901-91 30 (1 230) 


VJ 1 VJ/ 1 / \VJ 


ATA^AG 


^G^AG 




Non-codinQ rpQion(NCR) 


91 75_933i (1 57) 


9187-9277 (91) 

2 1 U/ J z- / / 


9131_923i (101) 

J \ ~J \ J A. ~J \ \ \ VJ 1 J 










tRNA-Mpt (M) 


9345-9407 (63) 


9288-9352 (65) 


9232-9294 (63) 








CAT 


tRNA-Gly (G) 


9407-9462 (56) 


9355-9410 (56) 


9296-9351 (56) 








TCC 


cox2 


9463-10158 (696) 


9411-10106 (696) 


9364-10050 (687) 


TO^AG 


TO^AG 


ATG^AG 




tRNA-His (H) 


10163-10218 (56) 


10113-10167 (55) 


10049-10105 (57) 








GTG 


rmi 


10211-11166 (956) 


10167-11121 (955) 


10103-11055 (953) 










nad3 


11167-11502 (336) 


11119-11454 (336) 


11053-11391 (339) 


TO^AA 


ATA^AA 


TO^AA 




nod5 


11502-13082 (1581) 


11451-13031 (1581) 


11394-12968 (1575) 


ATA^AA 


ATA^AG 


ATT/TAG 




tRNA-Ala (A) 


13085-13139 (55) 


13031-13085 (55) 


12968-13022 (55) 








TGC 


tRNA-Pro (P) 


13140-13197 (58) 


13094-13149 (56) 


13026-13082 (57) 








TGG 


tRNA-Val (V) 


13197-13252 (56) 


13149-13205 (57) 


13082-13137 (56) 








TAC 


nod6 


13253-13687 (435) 


13206-13640 (435) 


13138-13572 (435) 


m^AA 


TO^AA 


TTGfTAG 




nad4L 


13702-13941 (240) 


13656-13895 (240) 


13588-13827 (240) 


ATT/TAA 


GTT/TAA 


ATT/TAA 




tRNA-Trp (W) 


13919-13976 (58) 


13873-13930 (58) 


13805-13862 (58) 








TCA 
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Ascaridia galli;Ascaridia co!umbae;Ascaridia sp. (Ascaridida) 



GAl 
GA2 
GA3 
GA4 
GA5 
GA6 
GA7 
GA8 

GA9 
GA10 
GA11 
GA12 
GA13 
GA14 
GA15 
GA16 
GA17 

GA18 
GA19 
GA20 
GA21 
GA22 
GA23 
GA24 
GA25 



CO* I -C-N-Y-/K,</| -<///>f>-K-l .2-S I -/W2-l-R-Q-0-n-/77jS-S2-F-cTfb-l . I i-ox3-T-mui4-M-G-coA 2-H-rrnl ^naJ*-miJ*-A-P-V-muK>-ntuM.-\\ ! 



Heterorhabditis bacteriophora (Heterorhabditidae) 



vox I -\1-C-t ttv2-l l-/ /7»L-/M</3-Q-F-tT/b-L I -c av3-T-/irt</4-D-Ci-H t A/5-A-PAwwt/6-;wJ4L-W-E-m;S-S2-N-Y-/wt/l -<j//>6-K-L2-S 1 -mnfl-l-R 



Ancylostoma duodenale (Ancylostomatidae); Ascaris suum (Ascarididae) 32 species 



a I -C-M-D-G-c<>.v2-I l-r/7iL-/;<«/3-/w</5-A-P-V-/k/</6-/kK/4L-VV-E-r/7iS-S2-N-Y-/w < /l -<///>6-K-L2-S I -/k«/2-I-R-Q-F-c.» Ib-L I -mv3-T-/w</4 



Steinernema carpocapsae (Steinernematidae) 



:wv 1 -C-M-D-G-c K«2-l l-rrnL-muft-ikufi- A-PAwwJ6-N W4L-W-E-/rwS-S2-Y-//</</l -<///>6-K-L2-S I -w</<£-I-R-Q-F-n /b-L 1 -m\3-T-/wJ4 



Metastrongylus pudendotectus (Metastrongylidae) 



a I -( "-M-D-G-c ox2-\ I-r/7fL-/; < «/3-/;<«/5-A-I-P-V-/;<«/6-/;w</4L-W-E-m»S-S2-N-Y-/k«/l -<tf/>6-K-L2-S 1 -//</<£-R-Q-F-o/b-L I -COXi-T-nadA 



Cucullanus robustus (Cucullanidae) 



a I -C-M-D-G-i«a2-I I-miLww</3 W5 W6wW4LAV-A-E-miS-S2-YWl -<///*- K-L2-S I -/*/<£-I-R-Q-F-o tb-L I -m\3-T-W4-V-N-P 



Radopholus similis (Pratylenchidae) 



a I -P-RW6-W4L-V-K-S I -( , -L2-M-T-/w</4-FH>/b-D-;wJ2-Um3-N-E-(i-m>S-S2A'AV-/w</l -<i//>6-W5-Q-A-t <«2-I \-rmh-tuuti-L 1 



Chandlerella quiscali (Onchocercidae) 



a I -W-W6-R<^c v/b4 J ^y«3^^ 



Dirofllaria i mm it is (Onchocercidae); Loa loa (Onchocercidae) 7 species 



v I -W-w<i</6-R-Q.c;i/b-L 1 -coxh A-L2-N-M-K-/w</4L-m>S- Y-notil -F-<7//)6-l-G-mv2-H-m;L-/w < /3-C-S2-P-D- VW5-E-S I -/w</2-T-/w</4 



Onchocerca volvulus (Onchocercidae) 



v I -\VW6.R-Q-c.i/b-L 1 -<m3-K-A-L2-N-M-/W(/4L-r/7;S- Y-/k«/l ■F-<7//Hj-l-G-mv2-H-m;L-/w < /3-C-S2-l > -D-V-AW < /5-E-S I -thhtl-Y-mJA 



Heliconema longissimum (Physalopteridae) 



a I -W-W6-R-Q-C »/b-L I -t •Mv3-A-L2-K-V-K-/Mif/4L-r/7fS-YWI -F-<///>6-I-G-cym2-I I-r™L.-/j</</3-(-S2-P-D-M-/;</</5-E-S I -tuutl-J-mulA 



Wellcomia siamensis (Oxyuridae) 



a I -mu/l -<///>6-L2-S I -/w</24A'-R-Q-(V/7iLA1-m/</6A'AV-F^^ I -E-D-to.t2-I I-rmS-A-S2-N 



Enterobius vermicularis (Oxyuridae) 



vox I -iuui\ -<///*-L2-S I -mutl ■Y-R-Q-(^r/7>L-M-/KK/6A'-W-FH'>'/b^m3-P-TW4-(i W4L-K-/KftA3-w^M 1 -E-D-t <n2-I I-m/S- A-S2-I-N 



Dracunculus medinensis (Dracunculidae) 



vl -C-M-D-G-or>v2-H v77>L-n™flvw</5-A-P^viad6*fl^ 1 <»x3-T-0<KAt-W-S I -L-l-N-F-P-Y- VWI -K-L2-«//>6-T 



Strongyloides stercoralis (Strongyloididae) 



v I -nad4L-KS2-nad5-nad4-V/-R-S I -*-m-3-u//^-H-0-N-T-P-*«ii/2-V-l-/K»/3-D-iK«/ 1 - A-miS-M-C-rr/b-l . I -Cax2-E-\ 2-F-m;l .-/K/J6-G-T 



Agamermis sp. (Mermithidae) 



rl J-R^ Vw^U N^a«3^I^«rg^AJC^M</l -S 1 -h^wkkA- A-S2-wtS-C A S2 vM^ V.E-A-S2j>^i^^to64^4 J lvMrf5.W4?.Y.L2.P^.H.T 



Hexamermis aprotis (Mermithidae) 



or 1 -W-rtod 1 -D-l-R- /M<fc - cox3 -G.Y ^/p6- S 1 -S 1 -atD6-aiD<>-S 1 -E-S2-rrnS-Q-n(Kft-L2- nad\ - ftad4L - V-T-w<«/2-a1b-K-*r»L- W-F-H-M<< >r2-L 1 • A-C-wotfl -P 



Romanomermis culicivorax (Mermithidae) 



M K /■,;./< N..m./I-I-II-s:-Wm 



Romanomermis iyensari (Mermithidae) 



AUv/b-S2-WN\V^m-3-R-w</SV-C-AA^^ 



Romanomermis nielseni (Mermithidae) 



xUm^-R4-rmS-W4-K-M^/l-H-c < n-2-P-Y-M^/4i.-0->) < ; < /f>-1.2-» < <</2-S I -K-Sl -noiff -V-C- A-NfifpftWI -S2-c wb-H-SI -W-l. I -T-rwl.-M-D-G-F 



Strelkovimermis spiculatus (Mermithidae) 



iv 1 -D-W<t>y3-Y-c«-tf 2-F-R-<t//)6-L I -l-*<j(i&-K-/w</5-C-A-S I -m/L-mtS-Q-wrKfl-o fb-H AVw</4L-L2-E-P-(^/w</2vj<k/4-T-M-S2-W1 -N-C-A-S I -rrnL -m»S 



Thaumamermis cossrovei (Mermithidae) 



q>xl-C-A^Y>x3-WUv/b-w;</4<;^ 



Trichinella spiralis (Trichinellidae) 3 species 



iym1^w2-L2-E-w^I-K-W2A14-a?<<<^ 



Xiphinema americanum (Longidoridae) 



Bursaphelencbus xylophilus (Aphelenchoididae) 



Al-C-M-I>G<0x2-H-mfL-/w</3-//<^ 



Figure 2 Mitochondrial gene arrangement in three Ascaridia species (GAl) compared with those in other nematodes (GA2-GA25). 

The circular mt genomes were linearized at the 5' end of cox 7 gene for illustration purpose. Non-coding regions were not shown. 
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[Ascaridia sp. (GHL-2012)] respectively. A NCR at this 
location, i.e. between nad4 and trnM, is also present in 
the mt genomes of most other nematodes [27,28]. 

Gene rearrangement in the mitochondrial genomes of 
Ascaridia species 

A. galli, A. columbae and Ascaridia sp. (GHL-2012) have 
an identical gene arrangement in their mt genomes. The 
arrangement of mt genes in these three species, however, 
is substantially different from those of other nematodes 
(Figure 2). The mt genomes of the 62 species of nema- 
todes that have been sequenced to date showed 25 types 
of gene arrangements (GA1-GA25, Figure 2), of which 
GA3 is the most common type and has been found in 
32 species of nematodes. In comparison to the GA3 
type, at least three rearrangement events occurred in the 
Ascaridia species: trnM was translocated, and a block of 
11 genes (trnN, trnY, nadl, atp6, trnK, trnL 2 , trnS h nad2, 
trnly trnR and trnQ) swapped position with another block 
of 12 genes (trnG, cox2, trnH, rrnL, nad3, nad5, trnA, 
trnP, trnV, nad6, nad4L and trnW) (Figure 3). 

Phylogeny of nematodes inferred from mitochondrial 
genome sequences 

Monophylies of the order Rhabditida and the subclass 
Dorylaimia 

Of the 65 species of nematodes included in the phylo- 
genetic analyses in this study, 54 species were from the 
order Rhabditida of the class Chromadorea, and 11 spe- 
cies were from the subclass Dorylaimia of the class 
Enoplea. Both Rhabditida and Dorylaimia were mono- 
phyletic in all of the trees inferred with Bayesian, ML 
and MP methods from concatenated amino acid se- 
quences deduced from the sequences of the 12 mt 
protein-coding genes. The monophyly of the order 
Rhabditida was strongly supported with a posterior 
probability (PP) of 1 in Bayesian analysis (Figure 4), a 
bootstrapping frequency (BF) of 100% in ML analysis 
(Figure 5), and a BF of 94% in MP analysis (Figure 6). 
The monophyly of the subclass Dorylaimia was strongly 
supported in Bayesian and ML analyses (PP = 1, Figure 4; 
BF = 97, Figure 5), and was moderately supported in MP 
analysis (BF = 67%, Figure 6). 



Phylogenetic relationships within the order Rhabditida 

The 54 species of nematodes in the order Rhabditida 
included in this study were from three suborders: 
Rhabditina (24 species), Spirurina (26 species), and 
Tylenchina (4 species). The monophyly of the suborder 
Rhabditina was strongly supported in Bayesian and ML 
analyses (PP = 1, Figure 4; BF = 99%, Figure 5), and was 
moderately supported in MP analysis (BF = 79%, Figure 6). 
The Spirurina and the Tylenchina, however, were not 
monophyletic in all of the three phylogenetic analyses in 
this study. 

Of the 24 species from the suborder Rhabditina, 23 
species were from the infraorder Rhabditomorpha and 
one species from the infraorder Diplogasteromorpha. 
Species of these two infraorders were most closely 
related in all of the phylogenetic analyses in this study. 
The Rhabditomorpha, however, was paraphyletic with 
respect to the Diplogasteromorpha. Three species from 
the superfamily Rhabitoidea of the Rhabditomorpha 
were more closely related to Pristionchus pacificus 
(Diplogasteromorpha) than they were to the other 20 
species from the Rhabditomorpha. The close relation- 
ship between the species of the superfamily Rhabitoidea 
and P. pacificus was strongly supported in Bayesian ana- 
lysis (PP = 1, Figure 4), and was moderately supported 
in ML and MP analyses (BF = 83%, Figure 5; BF = 82%, 
Figure 6). In addition to the Rhabitoidea, three other 
superfamilies of the infraorder Rhabditomorpha were also 
represented in our analyses: Ancylostomatoidea (4 spe- 
cies), Strongyloidea (10 species), and Trichostrongyloidea 
(6 species). The Trichostrongyloidea was monophyletic 
with strong support in all of the three phylogenetic 
analyses in this study (PP = 1, Figure 4; BF = 88%, Figure 5; 
BF = 88%, Figure 6). The Rhabitoidea was monophyletic 
with strong support in Bayesian analysis (PP = 0.99) and 
weak support in MP analysis (BF = 52%, Figure 6), but was 
paraphyletic in ML analysis with weak support (BF = 51%, 
Figure 5). The Ancylostomatoidea and the Strongyloidea 
were paraphyletic in all of the three phylogenetic analyses 
in this study; there was strong support for the paraphylies 
of these two superfamilies in Bayesian and MP analyses 
(PP > 0.96, Figure 4; BF = 99%, Figure 5). 

The 26 species of the suborder Spirurina included in this 
study were from three infraorders: Ascaridomorpha (14 
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Figure 3 Rearrangement of mitochondrial genes in three Ascaridia species (pattern GA1) relative to the most common pattern of mt 
gene arrangement observed in nematodes (GA3). 
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Figure 4 Phylogenetic relationships among 65 species of nematodes inferred from Bayesian analysis of deduced amino acid 
sequences of 12 mitochondrial proteins. Taenia multiceps and T hydatigena (GenBank accession numbers FJ495086 and FJ518620) were used 
as the outgroup. Posterior probability (PP) values were indicated at nodes. 



species), Oxyuridomorpha (2 species) and Spiruromorpha 
(10 species). The Oxyuridomorpha and the Spiruromorpha 
were both monophyletic with strong support in all of the 
three phylogenetic analyses (PP = 1, Figure 4; BF = 100%, 
Figure 5; BF = 100%, Figure 6). The Ascaridomorpha, how- 
ever, was monophyletic only in ML analysis with moderate 
support (BF = 84%, Figure 5), and was paraphyletic in 
Bayesian analysis with strong support (PP = 0.99, Figure 4) 
and MP analysis with weak support (BF = 52%, Figure 6). In 



both Bayesian and MP analyses, the 10 species of the super- 
family Ascarodoidea of the infraorder Ascaridomorpha 
were more closely related to those of the Rhabditina and 
Steinernema carpocapsae (Tylenchina) than they were to 
the three Ascaridia species and Cucullanus robustus, which 
were also from the infraorder Ascaridomorpha. The three 
Ascaridia species, for which mt genomes were sequenced 
in this study, form a monoplyletic group with strong support 
in all of the three phylogenetic analyses (Figures 4, 5, 6). 
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Figure 5 Phylogenetic relationships among 65 species of nematodes inferred from maximum likelihood (ML) of deduced amino acid 
sequences of 12 mitochondrial proteins. Taenia multiceps and T. hydatigena (GenBank accession numbers FJ495086 and FJ518620) were used 
as the outgroup. Bootstrapping frequency (BF) values were indicated at nodes. 



Cucullanus robustus was sister to Ascaridomorpha + 
Rhabditina + Steinernema carpocapsae in Bayesian ana- 
lysis with strong support (PP = 1, Figure 4). Cucullanus 
robustus, however, was sister the Ascaridia species with 
moderate support in ML analysis (BF = 83%, Figure 5) and 
weak support in MP analysis (55%, Figure 6). 



The four species of the suborder Tylenchina included in 
this study were from two infraorders: Panagrolaimorpha 
(2 species), and Tylenchomorpha (2 species). Both of these 
infraorders were paraphyletic in all of the three analyses 
(PP>0.99, Figure 4; BF > 57%, Figure 5; BF>51%, 
Figure 6). 
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Figure 6 Phylogenetic relationships among 65 species of nematodes inferred from maximum parsimony (MP) of deduced amino acid 
sequences of 12 mitochondrial proteins. Toenio multiceps and T. hydotigeno (GenBank accession numbers FJ495086 and FJ518620) were used 
as the outgroup. Bootstrapping frequency (BF) values were indicated at nodes. 



Phylogenetic relationships within the subclass Dorylaimia 

The 11 species of the subclass Dorylaimia included in 
this study were from three orders: Mermithida (7 spe- 
cies), Dorylaimida (1 species), and Trichinellida (3 spe- 
cies). The Mermitida and the Trichinellida were both 
monophyletic with strong support in all of the three ana- 
lyses (PP = 1, Figure 4; BF = 100%, Figure 5; BF = 100%, 
Figure 6). Among the three orders, the Mermithida and 
the Dorylaimida were more closely related than either of 
them to the Trichocephalida, with strong support in 
Bayesian analysis (PP = 1, Figure 4) and moderate support 



in ML analysis (BF = 73%, Figure 5). The Dorylaimida was 
grouped with the Trichocephalida in MP analysis but the 
support to this grouping was weak (BF = 40%, Figure 6). 

Discussion 

Phylogenetic relationships among nematodes have been 
revised in the past two decades using sequences of 
nuclear small subunit (SSU) rRNA gene [12-14]. The 
current working hypothesis of Nematoda phylogeny 
incorporated evidence from both morphology, ecology 
and nuclear SSU rRNA gene sequence analyses [14,15]. 
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Recently, Kang et al. [23] tested the hypothesis of the 
Nematoda phylogeny with the mt genome sequences of 
25 species. Kang et al. [23] showed strong support to 
the monophyly of the class Chromadorea, represented 
by 16 species from the order Rhabditida; this study, 
however, could not establish the monophyly of the class 
Enoplea, represented by nine species from the subclass 
Dorylaimia. Park et al. [24] and Sultana et al. [25] ex- 
panded the taxon sampling and inferred phylogenies 
with mt genome sequences of 36 and 41 species of 
nematodes respectively. Park et al. [24] showed strong 
support for the monophylies of both Chromadorea and 
Enoplea. Sultana et al. [25], however, could not establish 
the monophyly of the class Enoplea in most of their 
analyses. A major difference between the mt genome 
phylogenies and the current working hypothesis of the 
Nematode phylogeny is on the monophyly of the suborder 
Spirurina, which includes Clade III nematodes proposed 
in Blaxter et al. [12]. In the current working hypothesis, 
the suborder Spirurina is monophyletic, which includes 
the infraorders Ascaridomorpha, Spiruromorpha and 
Oxyuridomorpha, whereas the infraorder Rhabdito- 
morpha is in the suborder Rhabditina [13,14]. Further, 
Ascaridomorpha is most closely related to Spiruro- 
morpha, whereas Oxyuridomorpha is sister to the group 
that includes Ascaridomorpha + Spiruromorpha. In the 
phylogenies inferred from mt genome sequences, 
however, Ascaridomorpha is most closely related to 
Rhabditomorpha; in most analyses, Oxyuridomorpha is 
sister to Ascaridomorpha + Rhabditomorpha, whereas 
Spiruromorpha is sister to the group that contains 
Ascaridomorpha, Rhabditomorpha and Oxyuridomorpha. 
Kang et al. [23], Park et al. [24] and Sultana et al. [25] 
contained much less taxa compared to previous phylogen- 
etic studies on nematodes with nuclear SSU rRNA gene 
sequences (e.g. >200 taxa in Meldal et al. [14]). The novel 
phylogenetic relationship inferred from mt genome se- 
quences, thus, needs to be tested further with more taxa 
from a wide range of nematode lineages. 

In the present study, we sequenced the complete mt 
genomes of three Ascaridia species from the family 
Ascaridiidae of the superfamily Heterakoidea, which was 
not represented in previous studies (Kang et al. [23]; 
Park et al. [24]; Sultana et al. [25]). Furthermore, we 
inferred the phylogenetic relationships among 65 nema- 
tode species, for which complete mt genomes have been 
sequenced to date. Our analyses support the division of 
nematodes into two classes, Chromadorea and Enoplea, 
represented by the order Rhabditida and the subclass 
Dorylaimia, respectively, in the present study. Both the 
Rhabditida and the Dorylaimia were monophyletic with 
strong support, regardless the methods of phylogenetic 
analysis used. Within the order Rhabditida, the suborder 
Rhabditina was monophyletic; the other two suborders, 



Spirurina and Tylenchina, however, were both para- 
phyletic in all of the three phylogenetic analyses in this 
study. Monophyly of the suborder Spirurina was well 
supported in nuclear SSU rRNA gene sequence analyses 
[12,14], but was rejected in mt genome sequence analyses 
in Kang et al. [23], Park et al. [24], Sultana et al. [25] and 
the present study. Phylogenetic analyses of mt genome 
sequences in all of these four studies support strongly a 
close relationship between Ascaridomorpha and Rhab- 
ditomorpha to the exclusion of Oxyuridomorpha and 
Spiruromorpha. Apparently, additional markers for phylo- 
genetic inference are required to resolve the controversy 
on the Spirurina between mt genome sequence phylogeny 
and nuclear SSU rRNA phylogeny. 

Gene arrangement in mt genomes provides a source of 
information for phylogenetic inference, independent 
from mt genome sequences [18]. It is noteworthy that 
10 species of Ascaridomorpha and 20 species of Rhab- 
ditomorpha included in our phylogenetic analyses share 
a common pattern of gene arrangement, GA3, in their mt 
genomes (Figure 2). GA3 is also present in Pristionchus 
pacificus, the only species from the infraorder Diplogas- 
teromorpha, which is closely related to the species from 
the superfamily Rhabitoidea of the Rhabditomorpha 
(Figures 4, 5, 6). Indeed, GA3 is the most common pattern 
of mt gene arrangement observed in nematodes to date; 
nevertheless, this pattern is present only in species of the 
Ascaridomorpha, the Rhabditomorpha and the Diplogas- 
teromorpha. The possibility that GA3 is ancestral to the 
phylum Nematoda can be rejected with confidence thanks 
to the discovery of mt gene arrangement pattern GA23 in 
three Trichinellida species [29,30]. GA23 shares obvious 
similarity with that of many arthropods and would resem- 
ble the ancestral pattern of mt gene arrangement of the 
phylum Nematoda more than any other patterns found in 
the nematodes. There is no evidence either that GA3 is 
ancestral to the class Chromadorea. Given the strong 
support to the close relationship between Ascaridomorpha 
and Rhabditomorpha + Diplogasteromorpha indicated by 
mt genome sequence analyses in the present study, 
the most parsimonious explanation is that GA3 is a shared 
derived character (i.e. synapomorphy) of the clade 
Ascaridomorpha + Rhabditomorpha + Diplogasteromorpha. 
The other six patterns, GA1, GA2, GA4, GA5, GA6 
and GA25, observed in either Ascaridomorpha and 
Rhabditomorpha or in species closely related to these 
two infraorders, can be derived from GA3 by only a few 
rearrangement events (Figure 2), e.g. GA1, observed in 
the three Ascaridia species, can be converted from GA3 
by three rearrangement events (see above). 

Monophyly of the suborder Tylenchina could not be 
established previously in nuclear SSU RNA gene se- 
quence analyses [14]. No species from this suborder 
were included in two previous analyses of mt genome 



Liu et al. BMC Genomics 2013, 14:414 
http://www.biomedcentral.com/1471 -21 64/1 4/41 4 



Page 10 of 13 



sequences [23,24]. The present study included four 
species from the Tylenchina and indicated the paraphyly 
of this suborder. The positions of the four species of the 
Tylenchina change depending on the methods of phy- 
logenetic inference used in this study. For instance, 
Steinernema carpocapsae is sister to the suborder 
Rhabditina in Bayesian analysis (Figure 4) and MP ana- 
lysis (Figure 6); it is, however, sister to Ascaridomorpha + 
Rhabditomorpha + Diplogasteromorpha in ML analysis 
(Figure 5). Further test with more species from this 
suborder apparently is needed in order to clarify the 
phylogenetic position of the suborder Tylenchina. Indeed, 
a very recent study by Sultana et al. [25] included six 
species from the Tylenchina and showed the paraphyly of 
this suborder with strong support. The paraphyly of 
the infraorder Rhabditomorpha with respect to the 
Diplogasteromorpha will also need further test with ex- 
panded taxon sampling as only one species, Pristionchus 
pacificus, from the Diplogasteromorpha was included in 
the present study. No species from the subclass Enoplia, 
which has two orders and eight suborders [13], was in- 
cluded in our analyses. For the subclass Chromadoria, 
all of the species included in this study were from the 
order Rhabditida whereas the other five orders, Plectida, 
Araeolaimida, Monhysterida, Desmodorida, and Chro- 
madorida, were not represented. Expanding taxon sam- 
pling from these lineages of nematodes is clearly the 
next step for phylogenetic studies of nematodes with mt 
genome sequences. 

Conclusions 

We sequenced the complete mt genomes of three 
Ascaridia species (family Ascaridiidae) that infest chick- 
ens, pigeons and parrots, respectively. These species have 
an identical arrangement of mt genes to each other but 
differ substantially from other nematodes. Phylogenetic 
analyses of the mt genome sequences of the Ascaridia 
species, in combination with 62 other nematode species, 
support the monophylies of the subclass Dorylaimia; the 
orders Rhabditida, Trichinellida and Mermithida; the 
suborder Rhabditina; and the infraorders Spiruromorpha 
and Oxyuridomorpha. Analyses of mt genome sequences, 
however, reject the monophylies of the suborders 
Spirurina and Tylenchina, and the infraorders Rhabdi- 
tomorpha, Panagrolaimomorpha and Tylenchomorpha. 
Monophyly of the infraorder Ascaridomorpha varies 
depending on the methods of phylogenetic analysis. The 
Ascaridomorpha was more closely related to the infra- 
orders Rhabditomorpha and Diplogasteromorpha (sub- 
order Rhabditina) than they were to the other two 
infraorders of the Spirurina: Oxyuridorpha and Spiruro- 
morpha. The closer relationship among Ascaridomorpha, 
Rhabditomorpha and Diplogasteromorpha was also sup- 
ported by a shared common pattern of mitochondrial 



gene arrangement. Analyses of mitochondrial genome se- 
quences and gene arrangement in the current study and 
several previous studies have provided novel insights into 
the phylogenetic relationships among several major line- 
ages of nematodes. Many lineages of nematodes, however, 
are underrepresented or not represented in these analyses. 
Expanding taxon sampling is necessary for future phylo- 
genetic studies of nematodes with mt genome sequences. 

Methods 

Ethics statement 

This study was approved by the Animal Ethics Committee of 
Lanzhou Veterinary Research Institute, Chinese Academy of 
Agricultural Sciences (Approval No. LVRIAEC201 1-009). 
The animals, from which the nematodes were collected post- 
mortem, were handled in strict accordance with good animal 
practices required by the Animal Ethics Procedures and 
Guidelines of the People's Republic of China. 

Sample collection and DNA extraction 

Adult Ascaridia galli were obtained from naturally 
infected chickens in Hunan province, China. Adult A. 
columbae and an undescribed Ascaridia sp. [hereafter 
Ascaridia sp. (GHL-2012)] were obtained from naturally 
infected pigeons and parrots in Guangdong province, 
China. These nematodes were washed in physiological 
saline, identified preliminarily to species based on host 
preference, morphological characters and predilection 
sites [31-33], fixed in 70% ethanol, and stored at -20°C. 

Total genomic DNA was isolated from individual nema- 
todes using sodium dodecyl sulphate/proteinase K treat- 
ment, followed by spin-column purification (Wizard® SV 
Genomic DNA Purification System, Promega). In order to 
independently verify the identity of the specimen, the 
region spanning ITS-1, 5.8S rRNA gene and ITS-2 was 
amplified from each individual nematode by PCR using 
previously reported primers [34] and sequenced directly. 

PCR amplification of the mitochondrial genomes of 
Ascaridia species 

We amplified fragments (400-500 bp) of coxl, nadl and 
nad4 genes by PCR using conserved primers (Table 2). 
The amplicons were sequenced from both directions 
using BigDye terminator v3.1 on ABI PRISM 3730 
platform. We then designed specific primers from these 
fragments for long PCR amplification (Table 3). We 
amplified the entire mt genome of individual specimens 
of A. galli, A. columbae and Ascaridia sp. (GHL-2012) 
by long PCR in three or four overlapping fragments, 
respectively. The three overlapping long-PCR fragments 
for A. galli were between coxl and nadl (-2.5 kb), be- 
tween nadl and nad4 (-7.5 kb), and between nad4 and 
coxl (-5.0 kb). The four overlapping long-PCR frag- 
ments for A. columbae and Ascaridia sp. (GHL-2012) 
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Table 2 Primers used to amplify short-PCR fragments 
from Ascaridia galli, Ascaridia columbae and Ascaridia sp 



Name of primer 


Sequence (5' to 30 


Reference 


For A. galli 






pcoxl 






AGCoxlF 


GAAGmGTAmGACTGGTAAGAA 


This study 


AGCoxIR 


CAGTGAGACCACCAATAGTAAACAA 


This study 


pnc/cfl 






JB11 


AGATOGTAAGGGGCCTAATA 


[36] 


JB12 


ACCACTAACTAATOACmC 


[36] 


pnad4 






AGNad4F 


CTOTOmAAl I I I I I ATGCTGCT 


This study 


AGNad4R 


AAGCGGCTAAAGCCTOGCATCACT 


This study 


For A. columbae 
and Ascaridia sp. 






pcox2 






ACCox2F 


mAAGmGTOTATOTOTGGm 


This study 


ACCox2R 


AAAATACACCAACCACAGGAAAACT 


This study 


pnad'\ 






JB11 


AGATOGTAAGGGGCCTAATA 


[36] 


JB12 


ACCACTAACTAATOACmC 


[36] 


pcox3 






ACCox3F 


TGGTA^CTGGAC I I I I I I I GAT 


This study 


ACCox3R 


CCAAACTACATCTACAAAATGCCAA 


This study 



were between coxl and nadl (-2.5 kb), between nadl 
and nad3 (~5 kb), between nad?> and cox2 (-2.5 kb) and 
between cox2 and coxl (-5.0 kb). Each long- PCR 
reaction (25 ul) contains 2.5 uL 2 mM MgCl 2 , 4.0 uL 
0.2 mM each of dNTPs, 2.5 uL 2.5 ul lOx rTaq buffer, 
0.25 uL 2.5 uM of each primer, 0.25 uL 1.25 U rTaq 
polymerase (Takara), and 2 ul of total genomic DNA. The 
PCR conditions were: 92°C for 2 min (initial denatu- 
ration), then 92°C for 10 sec (denaturation), 48-5 1°C for 
30 sec (annealing) and 60°C for 10 min (extension) for 
10 cycles, followed by 92°C for 10 sec, 48-5 1°C for 30 sec, 
and 60°C for 10 min for 20 cycles and a final extension at 
60°C for 10 min. PCR amplicons were column-purified 
and then sequenced using a primer walking strategy [35]. 

Sequence assembly and mitochondrial genome 
annotation 

Sequence reads of A. galli, A. columbae and Ascaridia 
sp. (GHL-2012) were assembled with ContigExpress 
program of the Vector NTI software package version 6.0 
(Invitrogen, Carlsbad, CA). The mt genome sequences 
of the Ascaridia species were aligned with those of other 
nematode species available in GenBank using Clustal X 
1.83 [37] to infer gene boundaries. The open-reading 
frames and codon usage profiles of protein-coding genes 
were analysed by the Open Reading Frame Finder 



Table 3 Primers used to amplify long-PCR fragments from 
Ascaridia galli, Ascaridia columbae and Ascaridia sp 



Name of primer Sequence (5' to 3') 



For A. galli 




AGC1N1F 


ATAGAAGmGTAmGACTGGTAAGAAGGAGGT 


AGC1N1R 


CACAATACCAGTAACCAAAGTAGCATAAACAG 


AGN1N4F 


TAAGTOTOAAGAAGGAGCAGGAGAGT 


AGN1N4R 


CAAAAATGGAAAAGAACACAAAGCAGCA 


AGN4C1F 


I I I I IAIGCIGU I IGIGI I C_ I I I ICCA 


AGN4C1R 


CCAAATAAAGTOCCAGCCACCTAAA 


For A. columbae 




ACC1N1F 


AGTOGACTGmATCCGCCmGA 


ACC1N1R 


AmCATAAGACACTCTCTGACCTC 


ACN1C3F 


GCCAGAGGTCAGAGAGTGTCTOTG 


ACN1C3R 


C^GCITCACTATACTCTA^GCCTGT 


ACC3C2F 


CAGGCAATAGAGTATAGTGAAGCAA 


ACC3C2R 


ATAGAAGGCACAGCCCAAGAATGAA 


ACC2C1F 


TAGTATGTGATGmGGGAATGC^ 


ACC2C1R 


C^ACACCAGTAGGCACAGCGAT 


For Ascaridia sp. 




ACC1N1F 


AGTOGACTGmATCCGCCmGA 


ACC1N1R 


AmCATAAGACACTCTCTGACCTC 


ASN1C3F 


GGTCAGAGGGmCTOTGAGATOCT 


ASN1C3R 


CAACCGAACAATCmATOCTCAACA 


ASC3C2F 


TGTOAGTAATAAAGATOTOGGTO 


ASC3C2R 


ACCAGGAATGTCACCATACTCATAACT 


ASC2C1F 


TGGAGGTOATAATCG 1 1 G 1 G 1 1 1 1 GC 


ASC2C1R 


CACAACAACAAAGGCTGAAAAACCATC 



(http://www.ncbi.nlm.nih.gov/gorf/gorf.html) using the 
invertebrate mitochondrial genetic code, and subse- 
quently compared with those of Ascaris suum [38] and 
Contracaecum rudolphii B [28]. Sequences of protein- 
coding genes were translated into amino acid sequences 
using the invertebrate mitochondrial genetic code in 
MEGA 5.0 [39]. Translation initiation and termination 
codons were identified by comparison with those of the 
nematodes reported previously [27,28]. The secondary 
structures of 22 tRNA genes were predicted using 
tRNAscan-SE [40] with manual adjustment [41]. Tan- 
dem repeats in the non-coding regions were found 
using Tandem Repeat Finder program (http://tandem. 
bu.edu/trf/trf.html) [42]. 

Phylogenetic analyses 

Considering the high degree of intraspecific variation in 
nucleotide sequences of mt genes of nematodes [43], we 
used the deduced amino acid sequences of mt proteins 
for phylogenetic analyses. Amino acid sequences were 
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deduced from the sequences of the 12 mt protein- 
coding genes common for all of the 65 nematode species 
and two outgroup species included in this study. The 
deduced amino acid sequences of each of the 12 mt pro- 
teins of these species were aligned individually and then 
were concatenated into a single alignment. Ambiguous 
sites and regions in the alignment were excluded using 
Gblocks (http://molevol.cmima.csic.es/castresana/Gblocks_ 
server.html) [44] with the default parameters. Two Taenia 
cestodes (T. multiceps and T. hydatigena, GenBank acce- 
ssion number FJ495086 and FJ5 18620, respectively) were 
used as the outgroup. 

Phylogenetic analyses were conducted using three 
methods: Bayesian, maximum likelihood (ML), and max- 
imum parsimony (MP). The MtArt model of amino acid 
evolution [45] was selected as the most suitable model 
of evolution by ProtTest 2.4 [46] based on the Akaike in- 
formation criterion (AIC). Bayesian analysis was with 
MrBayes 3.1.1 [47]. As MtArt model is not implemented 
in the current version of MrBayes, an alternative model, 
MtREV, was used in Bayesian analysis. Four independent 
Markov chains were run for 1,000,000 metropolis - 
coupled MCMC generations, sampling a tree every 100 
generations. The first 2,500 trees represented burn-in, 
and the remaining trees were used to calculate Bayesian 
posterior probabilities (PP). Maximum likelihood (ML) 
analysis was performed using PhyML 3.0 [48] with the 
MtArt model. 100 bootstrap replicates were run and 
bootstrapping frequencies (BF) were calculated. MP 
analysis was with PAUP* 4.0b 10 [49]; indels were treated 
as missing characters. A total of 1,000 random addition 
searches using tree bisection-reconnection (TBR) branch 
swapping were performed for each MP analysis. Boot- 
strapping frequency (BF) was calculated from 1,000 
bootstrap replicates with 10 random additions per repli- 
cate in PAUP. Phylograms were drawn using FigTree 
v.1.31 (http://tree.bio.ed.ac.uk/software/figtree/). 
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